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Abstract In this work we show a class of oscillating configurations for the evolution of the 
domain walls in Euclidean space. The solutions are obtained analytically. Phase transitions are 
achieved from the associated fluctuation determinant, by the decay rates of the false vacuum. 


Introduction. — The interest in studying nonlinear 
systems has gradually grown since the 1950s, due to the 
fact that nonlinearity is present in many different areas 
of science, including particle physics, plasma physics, cos¬ 
mology and field theory, among others. Nonlinear sys¬ 
tems, particularly those having solitonic excitations, play 
a prominent role in the modelling of several physical, 
chemical and biological systems as well. For instance, in 
the electric conductivity of organic materials, for which 
polarons and other polymer chain solitons provide con¬ 
ducting polymers [1]. Moreover, solitons are also related 
to electric conduction through DNA molecules [2]. In par¬ 
ticular, solitons play a fundamental role in field theory, 
encompassing applications from particle physics to con¬ 
densed matter. Non-abelian solitons are specially impor¬ 
tant in gauge theories of elementary particle physics [3]. 
Gravitating non-Abelian solitons were firstly discussed in 
the context of four-dimensional Einstein-Yang-Mills the¬ 
ory [4]. In the context of string theory, supersymmetric 
solitons further play an important role in the study of the 
non-perturbative sector and in understanding string du¬ 
alities. In condensed matter, the homotopy classification 
of finite energy configurations can be realized by similar 
forms as can be accomplished in field theory, when the 
solitary-wave solutions of SU(2) gauge theory are classi¬ 
fied likewise. 

When a model with a field dependent potential, which 
has two or more degenerate minima, is taken into account, 
different ground states arise at different portions of the 
space. For instance, the so called domain walls can be 
described [5,6], connecting different portions of the space 


where the field has different values for the potential degen¬ 
erate minima, even further in asymmetric scenarios [7] . In 
other words, the field configuration interpolates between 
the potential minima. Solitons can describe such field con¬ 
figurations presenting a localized and shape-invariant as¬ 
pect, and have a finite energy density [8]. The presence of 
those configurations is well understood in a wide class of 
models, that encompass monopoles, textures, strings and 
kinks as well [9]. 

In a seminal work by Coleman, a classical field theory 
was presented for the study of the false vacuum decay in 
theories involving asymmetrical potentials by analysing a 
kind of Xip 4 asymmetric model. In this theory, the rela¬ 
tive minimum corresponds to the false vacuum, whereas 
the absolute minimum corresponds to the true one. After 
that, Callan and Coleman presented the associated quan¬ 
tum corrections for the theory [10]. These preliminary re¬ 
sults paved an increasing interest in cosmology, providing 
models where the scalar field potential driving inflation 
has low energy minima, also known as false vacua. The 
absolute minimum of the energy density corresponds to 
the so called true vacuum state of the universe. 

An important application of this theory is the process 
that involves phase transitions in statistical mechanics. In 
this case, inside the false vacuum, it occurs the formation 
of a bubble of true vacuum which initiates the decay. From 
a quantum mechanical point of view, it corresponds to a 
tunnelling probability from a false vacuum to the true one. 
Moreover, the theory of vacuum decay has a considerable 
physical importance and can include the effects of grav¬ 
ity [11]. In fact, in a cosmological perspective, the early 
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universe had an extremely high energy density in a region 
of false vacuum. As it expands and cools down, it passes 
through a phase transition towards the true vacuum [12]. 

A wide class of problems involving phase transitions are 
supported by the </> 4 theory. However, due to the nonlin¬ 
earity of this theory it is very rare to analytically approach 
the problem. Recently an interesting model [15] was pro¬ 
posed, named asymmetrical double-quadratic model. This 
model is similar to the asymmetric A</> 4 characterized by 
the following potential 


R(^) = ^ 2 ^|0|-e^+i(e-l) 2 . (1) 


A similar model has been extensively studied in the lit¬ 
erature in several contexts. As an example, oscillons con¬ 
figurations [13], which are time-dependent and long living 
solutions, were obtained in the presence of the so-called 
signum-Gordon model [14]. 

On the other hand, some years ago the Asymmetrical 
Double-Quadratic (ADQ) model allowed the attainment 
of oscillating configurations in Minkowski space-time [16]. 
In fact, oscillating configurations were shown to be respon¬ 
sible for a delay in the transition from a false to a true vac¬ 
uum. Besides, a connection with phase transitions, which 
occur in ferromagnetic materials, was presented. In addi¬ 
tion, fermions bound states were proposed in a background 
of static solutions provided from ADQ model [17]. 

This paper is organized as follows: in Section 2 we 
present the asymmetrical model to be analysed and the 
classical field configurations. In Section 3 we obtain the 
analytical solutions for the scalar field. In Section 4 we cal¬ 
culate the respective decay rates. In Section 5 we present 
the conclusions. 


Classical field configurations. — The so called 
Generalized Asymmetrical Doubly Quadratic model 
(GADQM) was studied in [17], corresponding to a gen¬ 
eralized version of the ADQ model. The advantage of 
the GADQM is that the vacua can be chosen to repre¬ 
sent a slow-roll potential, yielding inflaton fields which 
are important in cosmological inflationary scenarios. The 
authors in [17] considered the problem of fermion bound- 
states and zero modes in the background of kinks of the 
GADQM. Motivated by that work, Brito, Correa, and Du- 
tra [18] showed an approach to construct nonlinear sys¬ 
tems with analytical multikink profile configurations. A 
prominent consequence of this work is the possibility of the 
resulting field configurations can be applied to the study of 
problems of condensed matter, cosmology, and braneworld 
scenarios. Thus, in order to employ this model, that more¬ 
over is exactly solvable, we will study the model whose 


potential is given by [17,18] 
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W) 

W) 
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— oo < <fi < —a , (2a) 


A2 
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— a < <fi < 0, 


(2b) 


A3 

2 



0 < 4> < a , 


(2c) 


A4 

2 



a < <j> < 00 , 


(2d) 


where Aj+i = Ai ^ and i = 1,2,3. Fig. 1 depicts 

this potential, presenting four asymmetric vacua located 
at ±3a/2, ±a/2. Moreover, we can see in Fig. 1 that 
the potential has three false vacua and a true vacuum. 
Thus, according to the pioneering work of Coleman [12], 
the false vacua can decay, where phase transitions play 
an important role in various phenomena. For example, it 
occurs in the nucleation process of statistical physics [19], 
in cosmological contexts [20], and in the Weinberg-Salam 
model [21] as well. 




Figure 1: Potential of the model (top) and corresponding field 
configuration (bottom) for v = 1, Ai = 1, a = 5, 61 = 1, 
62 = 0.5, &3 = 0.1 and 64 = 0. 


In this case it is possible to have a phase transition be¬ 
tween vacua. It is important to remark that a class of con¬ 
figurations in 1, 2 and 3 dimensions, in a similar potential, 
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was presented in refs. [16,18], proposing the existence of 
kink configurations in the presence of the ADQ model. At 
this point, it is important to point out that in refs. [16,18] 
the spacetime used by the authors is Minkowskian, and the 
field configurations correspond to the static case. Here we 
will work in a single scalar field theory in four-dimensional 
spacetime, which has the following Euclidean action 


where v = (N — l)/2. If the variable Z = p^f\j is taken 
into account, it yields the well known solution 

* j (Z) = a j K v (Z) + b j I v {Z)> ( 9 ) 

where K„(Z) and I V {Z) are modified Bessel functions of 
order v. Therefore, the complete solution has the form 





(3) 


where H(0) is given by the potential shown in Eqs.(2a)~ 
(2d). Note that the Euclidean theory is left invariant by 
0(4) transformations. Moreover, the Euclidean system is 
not similar as the Minkowskian system, since (a"g — x 2 ) 
for Minkowski coordinates is invariant, whereas the Eu- 

4 

clidean version is provided by (x M ) 2 . Hence, the two 

p=l 

systems are different, although the tunnelling processes in 
field theory are described by the field equations in Eu¬ 
clidean spacetime. Now, from the variation of the above 
action, the corresponding Euclidean equation of motion is 
provided by 


N 




(4) 


i—1 


where N represents the number of dimensions. More¬ 
over, as a consequence of the Euclidean theory, we can 
observe in the above equation an elliptic differential oper¬ 
ator, which opens new possibilities to find a different set of 
real non-singular solutions. For the model to be analyzed, 
Eq.(4) reads 


d 2 (j>j(p) , N d(j>j (p) 

~dp jMp) 


7 i j 1; ■ • ■ j 4 , 


(5) 

with p = 77 I 2 + t 2 , wherein x represents the three spa¬ 
tial coordinates and r the Euclidean time, A\ = 3a/2, 
A 2 = a/ 2 , A 3 = —a /2 and A 4 = —3a/2. 

Now, in order to solve Eq.(5), the transformations 
<f>j(p) = — Aj + ipj(p) in the fields <pj are performed. Thus, 
it is not difficult to conclude that 


d 2 Pj(p) N dipj(p) 

dp 2 p dp 


7^7) = °- 


( 6 ) 


Analytical solutions. — There are many possible 
ways to solve the above equation, each of them with ad¬ 
vantages and disadvantages. We will work with the one 
which makes possible to find real solutions. Thus, in order 
to solve the above equation, we choose 


<Pj{p)= P v ®j{p)- (7) 


77) — — Af + P v [ a jKv(p\/\i) + (10) 

Now we search for configurations where four different re¬ 
gions do exist, wherein the field is connected. At this 
point it is important to remark that the solutions in each 
region must be continuous through the p axis, what leads 

to 0i(pi) = 7(Pi)> 772 ) = 772) and <fo73) = (Mas)- 
On the other hand, we are looking for the case where 
0i (p = 0) = —3a/2 and 4>a{p —1 00 ) = 3a/2, which re¬ 
sults in 

77) =-y+P !/& i77a/A/), 0<p<pi (11) 

027) = ~ 2 d-p U [ a 2^ v {P'J\j)+b2lu{P\/\i)]i Pi < P < P2 

037) = 2 +^' l/ [ a 37ti/(p-\/7)+77(p\/A/)], P 2 < P < P3 

047) = y + P v a 4 ,K u (p y /\~ j ), p 3 < P < 00 . (12) 

To determine the constants aj and bj we use the condi¬ 
tion of continuity of the function 0j and of its derivative, 
in each point of the region p. Here, it is worth to remark 
that these conditions generate a set of six equations which 
allows the determination of the values of these constants. 
For instance, we show the profile of the field configuration 
in Fig. 2 for a given set of values of the parameters of the 
model. Moreover, it is necessary to impose that the first 


0(P) 



Figure 2: Field configuration for v = 1, Ai = 1, a = 5, b\ = 1, 
62 = 0.5, 63 = 0.1, 64 = 0, pi = 3, p 2 = 9 and ps = 20. 


Hence, the remaining function Qj(p) satisfies the Bessel 
equation 


2 d 2 77) 

P dp 2 


+ P 


d®j(p) 

dp 


(7 + A iP 2 )7(p) = 0, 


( 8 ) 


derivative of the field configuration is continuous at the 
intermediate points. Such a constraint produces transcen¬ 
dental equations which restrict the distances among the 
roots of the scalar field 0. 
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First-order phase transitions. — As argued in the 
seminal work by Coleman [12], a field theory with false 
vacuum allows to understand a large number of phenom¬ 
ena from condensate matter physics to cosmology. The 
key for this understanding comes from the decay prob¬ 
ability per unit volume, which is named decay rate. In 
particular, in a cosmological context, we can suppose that 
when the universe had been created, it was far from any 
vacuum state. As it has expanded and cooled down, it first 
evolved to a false vacuum instead of the true one. Thus, 
in such a scenario, when the cosmological time elapses, 
the universe should finally be settled in the true vacuum 
state. Therefore, in this section we calculate the decay 
rate of each phase transition that occurs from the false 
vacuum to the true vacuum. 

In fact, let us denote as usual the two stable static states 
by 4> = 4>f v (the false vacuum). Since <j>f v is unstable, its 
energy E((f>f v ) must present an imaginary part responsi¬ 
ble for its decay width, as Im E(<j>f v ) = —r(<//„)/2. In a 
weakly coupling framework, the width T (4>f v ) is exponen¬ 
tially small. The energy of the false vacuum reads a func¬ 
tional integral of type J T><pe~ s ^, where Se[4>] regards 
Eq.(3). Hence the width T, which corresponds to the de¬ 
cay rate, is determined by the contribution of a non-trivial 
saddle point, named the bounce 4>b, of the Euclidean ac¬ 
tion Se[(I>] [10-12]. 

It is important to remark that the decay rate represents 
the tunnelling probability from the false vacuum to the 
true one, per unit time per unit volume U, in such a way 
that in the semiclassical limit it assumes the form 


by the splitting rj(x) = Y2x Q’XPxix) ■ Hence the con¬ 
tribution of the false vacuum to the integral J T>(j)e~ s ^ 
reads Iq = JUx [da*] exp [—| Y2x ^ a i] > or equivalently, 
to I 0 = n A A- i/2 . The quantity ]~[ A A can be interpreted 
as the determinant of the operator — d^d^+V"((j)f v ), since 
it corresponds to the product of its eigenvalues. Thus 
Io = det[—+ V”(4>f v )\. Moreover, Eq. (14) can also 
present zero eigenvalues (modes). In fact, the center of the 
bounce may be located at any point in Euclidean space- 
time, leading to the existence of four zero modes around 
the bounce, determined by = S~ 1 ^ 2 d IJj (j)B{x), /i = 
0,1, 2, 3 [30]. By treating the zero modes in the same way 
as non-zero modes, then the corresponding integral would 
diverge [30], what imposes that the zero modes are dis¬ 
missed from the calculations in the previous formula for 
E c i , accordingly. 

Alternatively, the bounce action S E \<j>d\ is evinced in 
the right-hand side of Eq.(15). Hence, in general, the 
above mentioned pre-factor reads [23] 


SeI&ci} 2 
47T 2 ^ 


det[-d^d M + V"{<t>f v )\ 
det[—+ V"(4> c i)\ 


(15) 


Now, in order to evaluate the decay rate per unit volume 
and unit time T/U, firstly we need to find the classical 
solution (j) c i. Thus, from the Euclidean action (3), we have 
the following equation of motion 

+ t|) - V ' {M = °- (16) 


T/U ~ Kex-p(-S E {<f>ci]), (13) 

where the prefactor term K , which corresponds to quan¬ 
tum fluctuations with respect to the classical solution at 
the bounce field, is a coefficient that depends on the de¬ 
tailed form of V(<j>). Indeed, Coleman and Callan proved 
that the contribution of the bounce to the energy of the 
false vacuum reads [10] 


/SA[<M 2 \ SM,,.) 

detl-d^dfj, + V"(4>fv)} 

\ 

de^-cF 4 ^ + F"(0 c z)] 


In the above equation, the prime denotes the derivative 
with respect to the argument, cf>f v is the false vacuum, 
and (j) c i is a classical solution associated to the bounce 
c t>B■ Furthermore, the ringed determinant in Eq.(15) de¬ 
notes that the zero modes of the operator —d^d^ + V", 
corresponding to the translational invariance about the 
bounce location, are taken out. 

More precisely, the contribution of (f>f v is determined 
by a Gaussian integral over the perturbations about the 
bounce, namely 4>(x) = 4>f v + t]{x). Hence (3) yields the 
fluctuations S(rj) = f d 4 x [— d 2 + V"(<j>f v )] 77 . Ref. 

[30] introduces the eigenfunctions rjx as 

[-dl + V” (<j>f v )]rix = \vx , (14) 


Here, for the model here presented by Eqs. (2a-2d), 
< f> c i represents the classical solution for each region. Thus, 
using Eqs. (11-12) we can obtain the correct form of the 
Euclidean action, which is fundamental for the decay rate 
calculation. For simplicity, we can rewrite the Euclidean 
action in the following compact form 

4 

«c*] = E^ ( 1? ) 

3= 1 


where 


At) 


2tt u+1 

I>+1) 


rPj 


J2v+1 


' Pi-i 


1 

2 



2 

+ Vj(<i>j) 


dp, 


(18) 


with po = 0 and p^ = 00 . In fact, both the coefficients 
K and the exponential factor explicitly depend upon the 
action. Hence, for the model here analysed, the decay rate 
is given by the transition from the false vacuum, located at 
< j)^ = —3a/2, to the true vacuum, at (/>$ = 3a/2. Note 

that before reaching the point <pn), the decay process rolls 
through the other two metastable potentials, located at 
< f> ^ = —a/2 and (jy^ = a/2. Hence, the decay rate from 
the false vacuum to the true vacuum becomes 

3 

r/u = EC^), (19) 

j =1 
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with the following definition 


ri 

Uj 


1 

(2^F 


E n(n) q (n 
d E 


\n,m= 1 


x n^pl^ 51 ^). 

9=1 


det(Ad(j)) 

det(M^) 


( 20 ) 


where we used the relation in Eq. (18). Moreover, we 
employ the fluctuation operator ■M(j) in the background 
of the classical solution and its counterpart in the 

false vacuum, which are given by 

m u) = -d»d, + v;' + M 0+1) \ (2i) 

M° U) = -d^ + Vj'^W), ( 22 ) 

where <pG) —► 4$, for a = 1, 2,3,4. Now, it is convenient 
to use the fact that we have the 0(4) spherical symmetry, 
in such a way that the operators -M(j) and Ad^ can be 
decomposed with respect to 0(4) angular momenta. Thus, 
we can separate the operators into partial waves, which 
can be written as radial operators in the following form 


M U,1) = 

d 2 

N d 

dp 2 

p dp 

M U = 

d 2 

N d 

dp 2 

P dp 


A_^V + v i+1 (p), (23) 
Z(/ 4 2) +V J (p), (24) 


where Vj(p) = V-'(<f>^). For simplicity, let us define 


_ det(Ad W) ) 
^-det(Ad“. 0 )- 


(25) 


In this case, G(j^ can be computed by the so-called 
Gelfand-Yaglom (G-Y) method [22], which efficiently and 
allows to find the determinant of an ordinary differen¬ 
tial operator without necessity of compute its eigenvalues. 
The G-Y method states that for the radial operators (23) 
and (24), Gjyq can be easily computed in the following 
form 


G U,i) ~ 


lim 

p—>00 


*u.o(p) 


P + 1) 2 


(26) 


where ^(j,j)(p) and ’F” j .^(p) are regular solutions of the 
equations 


Here, it is important to highlight that there are three im¬ 
portant types of eigenvalue associated to the fluctuation 
operator. The first is the 1 = 0 sector, which has a nega¬ 
tive eigenvalue mode of the fluctuation operator. In this 
case, it is responsible for the instability of the configura¬ 
tion and leads to decay. The second is the l = 1 sector, 
where there is a four-fold degenerate zero eigenvalue (zero 
modes) of the fluctuation operator. From a physical view¬ 
point, these four zero modes are equivalent to the Gold- 
stone modes, that exist due to breaking of translational 


invariance. Finally, the third type consists of the l > 2 
sectors, with positive eigenvalues. However, we emphasize 
that in our analysis the zero modes are removed from the 
fluctuation determinant. 

Our next step, in order to compute the decay rate, is 
to compute the explicit form of the Euclidean action. To 
accomplish it, it is necessary to determine the integral (18) 
for each region. Let us then write p = (/V A,-, in such a 
way that Eq. (10) reads 

0t(C) = — A? -(- C [flj-^(C) + bjI„{Q], (28) 

where dj = aj/y/Xj, and bj = bj/y/Xj. In this case, Eq. 
(18) can be rewritten as 


qti) _ ,i 

~ a C?>) 


rG/X\ 


' Cj—i 


x 


2i/+l 


1 f d(j)j 


- 2 {~f)+ v AMO] 


with the following redefinitions 
2tt v+1 


I (i+) 


X v aT{v + 1) : 


VjlMO] y (</b(C) + Gj) 2 + bj 


d(, 

(29) 

(30) 

(31) 


where Gj corresponds to the value in each vacuum. In ad¬ 
dition, the derivative dcpj/dC, is crucial in our calculation, 
yielding 

d ^=C[-XK,-iX) + XX-XO\- (32) 

Other important results to find the complete Euclidean 
action form are the integrals 

J dCC m K v (Q = 2- (l/+2) 7rC m_l/+1 csc(7r^) 

- C 2l T 


4 r 


TO — V +1 


TO V V + 1 




1=1 


X}, 


e 


(33) 


and 


J dCC m IX0 = 2-G+^)( m -’'+ 1 T ^ w +^ + 1 ^ 


xUrz 


{Q 2 },{U 2 ,T 2 },^ 


(34) 


where 


Qi s m + (- |, ) l + l | Rl = 1 + ( _+, Ti =Qi + l. 


Furthermore, the function 'Hrz [{Qi} , {7\b, 71} , C 2 /4] is 
the regularized generalized hypergeometric function and 
T{X) is the gamma function. Thus, by using of the above 
relations and the usual proceeding presented in [23,24], we 
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can obtain, after straightforward calculations, the value of 
the fluctuation determinant (20) in each region, and as a 
consequence the decay rates. In this case, the transition 
rates can be determined from a straightforward example. 
For 1 / = 1, Ai = 1, a = 5, &i = 1, &2 = 0.5, b 3 = 0.1, 
64 = 0 , pi = 10 ~ 3 , P 2 = 10~ 2 and p 3 = 10 _1 , we find that 

T 1 /V 1 ~ 10 - 25 , (35) 

Y 2 /V 2 ~ 10- 14 , (36) 

r 3 /F 3 - 10” 5 . (37) 

At this point, it is important to remark that here the ad¬ 
vantage is that in each region the model is exactly solvable. 
Hence the decay rate was obtained in an analytical form. 
Another important consequence is that the approach used 
can be applied to study the phase transitions in models 
where we have piecewise potentials. 

Concluding Remarks and Outlook. We investi¬ 
gated a class of oscillating configurations for the evolution 
of the domain walls in Euclidean space. We found in an 
analytical form the configurations of the field and the de¬ 
cay rate of the false vacuum. 

Based on cosmological numerical values we can further 
incorporate a better description of the phase transitions, 
as follows. In the cosmological context, the Universe is 
known to have different dominant dynamical components 
since its origin. From radiation and relativistic matter at 
early time, to late time dark energy component, passing 
through the matter-dominated scenario. In further work 
we intend to physically interpret the vacuum decays of our 
model as the transition of different dynamical eras of the 
universe, i.e., the transition from a radiation to a matter- 
dominated scenario will be analysed merely as a transition 
from a false to a true vacuum and the same will happen 
for the other dynamical transitions of the universe. To 
accomplish it, we will need the values for |.t| and r for 
each transition era of the Universe. 

From Fig. 2, it is possible to perceive our argumenta¬ 
tion. One can interpret each of the regions J, II, III 
and IV as the different dynamical eras of the Universe. 
While regions II , III and IV would stand for the radia¬ 
tion, matter and dark energy-dominated eras, the region I 
could account for the inflationary era (check, for instance, 
[25-28]). Note that the scalar field evolution also makes 
one able to predict the universe fate, through a deep anal¬ 
ysis on region IV. Such an analysis may, for example, 
predicts the universe to expand forever, yielding a Big- 
Freeze model [29]. 

The transition rates obtained offer a further interpreta¬ 
tion for the dynamical evolution of the Universe. In fact, 
even if the state of the early Universe was cold enough not 
to provide a thermal transition to the true vacuum state, a 
quantum decay from the false vacuum to the true vacuum 
may still be possible through the barrier mechanism. 
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